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A blow-up criterion of strong solutions to the 2D 
compressible magnetohydrodynamic equations 


Teng WANG* 


Abstract 

This paper establishes a blow-up criterion of strong solutions to the two-dimensional 
compressible magnetohydrodynamic (MHD) flows. The criterion depends on the 
density, but is independent of the velocity and the magnetic field. More precisely, 
once the strong solutions blow up, the A°°-norm for the density tends to infinity. 

In particular, the vacuum in the solutions is allowed. 

Keywords: compressible magnetohydrodynamic equations; blow-up criterion; Cauchy 
problem; vacuum. 


1 Introduction 

We consider the system for the two-dimensional viscous, compressible magnetohy¬ 
drodynamic (MHD) flows in the Eulerian coordinates as follows, 

Pt + di v(pu) = 0, 

< (pu)t + di v(pu (g) u) + V-P(p) = pAu + (p + A)V(divu) + (V x H ) x H, (1.1) 
H t — \7 x (u x H) = uAH, dwH = 0, 

where p = p(x, t ), u = (u 1 ,u 2 )(x, t ), H = (H 1 , H 2 )(x, t ), and 

P(p) = Rp^(R> 0,7 > 1 ) ( 1 - 2 ) 

are the fluid density, velocity, magnetic field and pressure, respectively. R is a positive 
constant. Without loss of generality, we assume that R = 1. The constant viscosity 
coefficients p and A satisfy the physical restrictions: 

p > 0, p + A > 0. (1.3) 

The constant u > 0 is the resistivity coefficient which is inversely proportional to the 
electrical conductivity constant and acts as the magnetic diffusivity of magnetic fields. 
We consider the Cauchy problem of (11.111 with the initial data 

(p,u,H)\ t=0 := (p 0 ,uo,H 0 ), x € M 2 , (1.4) 

and the boundary condition at the far fields 

(p, u, H) —»• (0,0,0), as \x\ —y oo. (1.5) 
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There have been large literature on the compressible MHD system m by many 
physicists and mathematicians due to its physical importance, complexity, rich phe¬ 
nomena and mathematical challenges, see [2ll5ll8l [T0l[TTl[TTl[T8llTi^ lM ll25ll29]l3Tll32] and 
the references therein. When the initial density is uniformly positive, the local exis¬ 
tence of strong solutions to the three-dimensional compressible MHD was proved by 
Vol’pert-Khudiaev m- Then, Kawashima [18] firstly obtained the global existence 
when the initial data are close to a non-vacuum equilibrium in H'^-no rm. In partic¬ 
ular, the theory requires that the solutions have small oscillations around a uniform 
non-vacuum state so that the density is strictly away from the vacuum. 

For general large initial data, the global well-posedness of classical solutions to the 
compressible MHD system remains completely open. One of the main difficulties is 
that even the initial density is absence of vacuum, one could not know whether the 
vacuum states may occur or not within finite time. The system (11.11) may degenerate 
in the presence of vacuum, which produces new difficulty in mathematical analysis. 
Therefore, it is interesting to study the solutions with vacuum for the compressible 
MHD system. In fact, the local well-posedness of strong solutions of the system m 
in three dimensions was established by Fan-Yu [5] when the initial density may contain 
vacuum. Recently, Li-Xu-Zhang [2T] proved the global existence of classical solutions 
with vacuum as far field condition to (ll.ll) - (|1.5l) in M 3 for regular initial data with 
small energy but possibly large oscillations. However, the two dimensional case is 
quiet different from the three dimensional case when the far field condition is vacuum. 
Precisely speaking, the difference between 2D and 3D is that, if a function u satisfies 
Vu E L 2 (R 2 ), it is impossible to imply that u E L P (R 2 ), for any p > 1, while if 
Vu E L 2 (R 3 ), then u E L 6 (R 3 ). More recently, Lv-Huang |23] proved the local existence 
of strong solutions for the Cauchy problem of the two-dimensional compressible MHD 
equations with vacuum as far field density by weighted energy estimate and Lv et 
al. [25] generalized the previous work of Li-Xu-Zhang [21] on two-dimensional case. 

Until now, all the global existence of solutions for the compressible MHD equations 
were obtained with some ’’smallness” assumptions on the initial data. However, there 
still remains a longstanding open problem: whether the strong (or classical) solutions 
to the compressible MHD system m can exist globally or not? Thus, it is important 
to study the mechanism of blow up and structure of possible singularities of strong 
solutions to the compressible MHD system (O). Along this direction, Xu-Zhang in 
[32] established the following Serrin’s type criterion to the three-dimensional isentropic 
compressible MHD equations: 



where T* E (0, oo) is the maximal time of existence for strong (or classical) solutions 
(p,u), L r w denotes the weak L r -space and r, s satisfies 

2 3 

—I— <1, 3 < r < oo. (1-7) 

S T 

Then Huang-Li in [16] extended the result of [32] to the non-isentropic case. 

Although the Serrin type criterion for the Cauchy problem of three-dimensional MHD 
flows has been well established by Xu-Zhang in [32]. However, the two-dimensional case 
with the initial density containing vacuum becomes more difficult since the analysis 
of [32] for 3D case depends crucially on the L 6 -bound on the velocity, while in 2D case, 
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the velocity may not belong to L P (R 2 ) for any p > 1. The main aim of this paper is to 
establish a blow-up criterion for the two-dimensional compressible MHD system. 

Before stating the main results, we first explain the notations and conventions used 
throughout this paper. For R > 0, set 

Br = {x € R 2 | |x| < i?} , f fdx = I fdx. 

J J K 2 

Moreover, for 1 < r < oo and k > 1, the standard homogeneous and inhomogeneous 
Sobolev spaces are defined as follows: 

f L r = L r (R 2 ), D k ’ r = D k ’ r (R 2 ) = (u € L 1 1 oc (M 2 )|V fe u € L r (R 2 )}, 

I D 1 = D 1 ’ 2 , W k ’ r = VF fc,r (R 2 ), H k = W k ’ 2 . 


Next, we give the definition of strong solution to CEO as follows: 

Definition 1.1 (strong solution) ( p,u,H ) is called a strong solution to (jl.ll) in R 2 x 
(0, T), if for some qo > 2 and a > 1, 

'p> 0, px a eC([0,T]-L 1 nH 1 nW 1 ’i°), p t £C([0,T];L 2 nl«), 

< (i t ,F)eC([0,T];D 1 nD 2 )nL 2 (0,T ; D 2 -»), H € C([0, T];H 2 ), (1.8) 

. (p 1/2 u t , H t ) € L°°(0, T; L 2 ), (u t ,H t ) € L 2 { 0, T; D 1 ), 

and (p,u,H) satisfies (11.11) a.e. in R 2 x (0, T), where 

x = (e T \x\ 2 ) 1 ^ 2 log 2 (e + |x| 2 ). (1.9) 

Without loss of generality, assume that the initial density po satisfies 

[ p 0 dx = 1, (1.10) 

Jr 2 

which implies that there exists a positive constant Nq such that 



Our main result can be stated as follows: 


1 

2 


( 1 . 11 ) 


Theorem 1.1 Let fl = R 2 . In addition to (11.101) and (11.111) . suppose that the initial 
data (p 0 ,u 0 ,H 0 ) satisfy, for any given numbers a > 1, q > 2, 

/ Po > 0, x a p 0 € L 1 0 H 1 0 W 1 ’ 9 , poul + p^eL 1 , 

{ ( uo,H 0 ) € D 1 0 D 2 , x a / 2 H 0 € L 2 , x a / 2 VH 0 € L 2 , 1 ' j 

and the compatibility condition 

- pAu 0 - (p + A)Vdivrt 0 + VP(p 0 ) - (V x H 0 ) x H 0 = p)j 2 g, (1.13) 

holds for some g € L 2 . 

Let ( p,u,H) be a strong solution to the Cauchy problem (11.11) . (11.51) and (11.41) . sat¬ 
isfying (US]). If T* < oo is the maximal time of existence, then 

linr ||p||l°c (0) t ; loo) = oo. (1.14) 
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Remark 1.1 Theorem 11.11 means the blow-up criterion (11.141) is independent of the 
velocity u and the magnetic field H. 

We now make some comments on the analysis of this paper. As previously mentioned, 
the L 6 -norm of velocity u plays crucial role in 3D MHD system, while for 2D, u may 
not belong to L P (R 2 ), for any p > 1. The key observation of this paper is that, if 
we restriction the initial data in a smaller space, i.e. x a po £ T 1 (M 2 ) for some positive 
constant a > 1 (see (|1. 121) 1. we can show that ux~ v belongs to L P0 (R 2 ), for some 
positive constant po > 1, and p £ (0,1] (see (13.521) 1. here x = (e + lx) 2 ) 1 / 2 log 2 (e + |x| 2 ). 
For this , we need to manipulate the weighted energy estimates throughout the Section 
[3] below. 

Before finishing the introduction, we recall some related works. The global existence 
of weak solutions to the system (11.11) was established by Hu-Wang [TOlITT] . If there is no 
electromagnetic field effect, (11.11) turns to be the compressible Navier-Stokes equations. 
For the blow-up criterion of the compressible Navier-Stokes equations, we refer to H2 
and the references therein. 

The rest of the paper is organized as follows: In the next section, we collect some 
elementary facts and inequalities for the blow-up analysis. The main result, Theorem 
11.11 is proved in Section [3l 

2 Preliminaries 

In this section, we will recall some known facts and elementary inequalities which 
will be used frequently later. 

The local existence of strong solutions when the initial density may not be positive 
can be proved in a similar way as in [19] (cf. [24]). 

Lemma 2.1 Assume that the initial data (pq,uq,Hq) satisfy (11.121) . Then there exist 
a small time T\ > 0 and a unique strong solution ( p , u, H ) in the sense of Definition 
11.11 to the Cauchy problem (11.11) - (11.41) in R 2 x (0, T). 

Next, the following well-known Gagliardo-Nirenberg inequality (see 2fij ) will be used 
later. 

Lemma 2.2 (Gagliardo-Nirenberg) For p £ [2,oo),(7 £ (l,oo), and r £ (2,oo), 
there exists some generic constant C > 0 which may depend on p , q , and r such that 
for f £ R 1 (R 2 ) and g € L«(. R 2 ) n D^R 2 ), we have 

C||/||2 2(R2) ||V/||^ ( 2 R2) , (2.1) 

||„|| < r \\ n \\<l{r-2)/{2r+q{r-2))u X7 n2r/(2r+q(r-2)) /o o', 

II#IIc(r 2 ) — ^ II v 51Iz/-(r 2 ) ‘ 

The following weighted L p bounds for elements of the Hilbert space D 1 (R 2 ) can be 
found in [22], Theorem B. 1]. 

Lemma 2.3 For m € [2,oo) and 8 £ (1 + m/ 2,oo), there exists a positive constant C 
such that we have for all v £ D 1,2 (R 2 ), 

1 772 \ 1/m 

(log(e + |x| 2 )) 9 dxj < C\\v\\ L 2 ( Bl ) + C'||Vu|| L 2 ( ]R 2 ). (2.3) 



4 













The following lemma was deduced in |2L)j . we only state it here without proof. 
Lemma 2.4 For x as in (II.9j) . suppose that p E L°°(R 2 ) is a function such that 

0 < p < Mi, M 2 < / pdx, px a € L\R 2 ), (2.4) 

J B Nt 

for IV* > 1 and positive constants Mi, M 2 , and a. Then, for r € [2, 00 ), there exists a 
positive constant C depending only on Mi, M 2 , a, and r such that 

(X 2 P ^ rdx ) ~ CN * + IIp^IIl 1 !® 2 )) (ll/^llz 2 (R 2 ) + l|Vu|| L 2 (R2) ) , (2.5) 

for each v E {u E L> 1 (M 2 )| p x l 2 v E L 2 (M 2 )} . 

Next, for V* 1 = (— 82 , <9i), denoting the material derivative of / by / = ft+u-V f. We 
state some elementary estimates which follow from (12.11) and the standard L p -estimate 
for the following elliptic system derived from the momentum equations in © 2 : 

A F = div (^pii — di v(H 0 H — -|i7| 2 l2)^ , (2-6) 

and 

/jAw = V x • (^pii - di v(H 0 H - ^|iif I 2 )^ , (2.7) 

where 

F = (2p + A)divu — P(p), uj = d\u 2 — d^u 1 ■ (2.8) 

The symbol 0 denotes the Kronecker tensor product, i.e. H 0 H = {H l W) 2 x 2 and 
” 12 ” denotes the 2x2 unit matrix. 

Lemma 2.5 Let ( p,u,H) be a smooth solution of M.l}) . Then for p > 2 there exists a 
positive constant C depending only on p,p, and A such that 

l|VF|| iP ( K 2) + ||Vu;||lp(r 2) < Cdlpullip^) + |||77||ViL|||x, P ( R 2)), (2.9) 

II-^IIlp(k 2 ) + 11 ^ 11 ip (r 2 ) A C (Hp“IIl 2 (r 2 ) + |||tf||VM||| L 2 (R 2 )) 

' (||Vu|| L 2( R 2) + ||T > ||l2( R 2)) ^ , (2.10) 

l|Vu|| LP(R 2) < C (||ptt|| L 2( R 2) + |||i7||Vi7||| L 2( R 2)) / p 

■ (I|V w IIl 2 (r 2 ) + II^’IIl 2 (r 2 )) + C\\P\\lp(m2)- (2-11) 

Proof: On the one hand, by the standard L p -estimate for the elliptic systems (see im 
(12.61) (12.7p yield (12.91) directly, which, together with (12.11) . (|2.6I) and (12.71) . gives (12.101) . 
On the other hand, since —An = —Vdivu — we have 

Vu = —V(—A) _1 Vdivu - V(-A)- 1 V- l cj. (2.12) 

Thus applying the standard L p -estimate to (12.121) shows 
HVullipptf) < C'(p)(||divu|| LP ( R 2 ) + |MIlp(r 2 )) 

< C'(p)||T||lp(r 2 ) + C(p)\Mlp(r2 ) + C(p)II-F1Ilp(r 2 )> 
which, along with (12.10ft . gives (12.111) . Then, the proof of Lemma 12.51 is completed. 

Finally, the following Beale-Kato-Majda type inequality, which was proved in m, 
will be used later to estimate ||Vu|| L°° and ||Vp|| L 2 ni9 (g > 2). 
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Lemma 2.6 For 2 < q < oo, there is a constant C(q ) such that the following estimate 
holds for all Vu € L 2 (M 2 ) n D 1 ' q (R 2 ), 

||Vu|| L =o (M 2 ) < C (||divtt|| L oo( R 2 ) + IM|l°°( R 2 )) log(e + ||V 2 u|| L 9 ( R 2 )) 

+ f 1 |V«|| L 2 ( R 2 ) + c. 


3 Proof of Theorem 11.11 


Let (p,u,H) be a strong solution to the problem (ll.ll) - (|1.4D as describe in Theorem 
11.11 Suppose that (11.141) were false, that is, 


I|PIU~(0,T;L~) = M 0 < OO. 


(3.1) 


First, the standard energy estimate yields 


sup 
o <t<T . 


1 . l2 P \H \ 2 

-p\ur -\ -h J — L 

2' 1 1 7-1 2 


^ dx + II (p\Vu\ 2 + v\VH\ 2 )dxdt<C, (3.2) 


for 0 < T <T*. Throughout this paper, several positive generic constants are denoted 
by C and C % (i = 1,2) depending only on Mo, p,, A, v. T*, q , a and the initial data. 
The following lemma is based on (13.21) . 


Lemma 3.1 It holds that for q £ [2,oo) and 0 < T < T*, 
\\H\\l°°(0,T;Li) + / 


\H\ q ~ 2 \S7H\ 2 dxdt < C. 


(3.3) 


Proof: We prove (13.31) as in He-Xin E3- Multiplying (II. IL by q\H\ q 2 H and inte¬ 
grating the resulting equation over M 2 yield that 


\H\ q dx + v J (q\H\ q ~ 2 \VH\ 2 + q(q-2)\H\ q ~ 2 \V\H\\ 2 )dx 
H ■ S7u ■ H\H\ q ~ 2 dx — (q — 1) J divu\H\ q dx 


= q 


< C||Vu|| L 2|||M|'?/ 2 || 2 4 < C||Vu|| L 2|||M|^ 2 || L 2||V|M|«/ 2 || i 2 
<<5||V|M| 9 / 2 ||| 2 +C , (<5)||Vu|| 2 r2 |||M| 


(3.4) 


Choosing S suitable small in (13.41) . we obtain (13.31) directly after using Gronwall’s in¬ 
equality and (13.21) . Thus the proof of Lemma l3Jl is completed. □ 


Next, we give the key estimate on Vu and VF in the following lemma. 
Lemma 3.2 Under the condition (13.11) . it holds that for 0 < T < T*, 



(3.5) 

(3.6) 
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Proof: First, multiplying the momentum equation O 2 by u and integrating the 
resulting equation over M 2 gives 


I p\u\ 2 dx = - J u- VPdx + fi J u - Ixudx + (p + \) J u- Vdivu dx 

5 

- ]- f u ■ V\H\ 2 dx + f u-H- VHdx = J2 Ii, 

•' •' i =1 


where we have used the fact that 

(V x H)x H = di v(H ® H) - ^S7\H\ 2 = H ■ VP — ^V|P| 2 . 
Since P satisfies 

P t + di v(uP) + (7 — l)Pdivtt = 0, 
integrating by parts yields that 


— u ■ VPdx = 


/ i 

< ( / di vuPdx^J +C||V 


U\\ 2 L 2. 


Integrating by parts also leads to 


p, ii ■ A udx = — 


(||Vu|| L 2 )* - / 2 J d i u j di(u-'Vu j )da 


(3.7) 


(3.8) 


J (^(divu)tP — (u ■ Vu) ■ VP^jdx 

di vuPdx^J + J ^(7 — l)P(divii ) 2 + PdiU ■ Vu^jdx ( 3 . 9 ) 


<-$(l|v«ll! 2 ) t + c||Vu||f 3 


(3.10) 


and that 


(/r + A) f u- Vdivu dx = — — ^ — (||divu|| 2 2 ) t — (/r + A) f divw div(rt • Vu)dx 


< -^r- (ll divM H h) t + C\\Vu\\ 3 L3 . 


(3.11) 
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Using (I1.1IU and (13.21) . we get 


h = 


i j \H\ 2 divutdx + i y |i/| 2 div(u • X7u)dx 
(/ ■* diwucfa^ + ^ J u ■ V\H\ 2 divudx + ^ J \H \ 2 div(u ■ Vu)dx 
— J (H ■ Vu + I'Aff — //divu) • i^divzidx 

^ divudx^ ~ J l^1 2 (divu) 2 dx + ^ J \H\ 2 diU ■ Vu l dx 


< 


< 


— J (H ■ Vu + z^Ai? — i/divu) • i7divrzdx 
J ^ di vudx^J +C(e) J \H\ 2 \Vu\ 2 dx + e\\V 2 H\\ 2 l2 
<1 f ^^divudx^j +C(e)\\Vu\\ 2 L 3 \\H\\% 3 \\V 2 H \\% 3 + e\\V 2 H \\ 2 L 2 

J —^-di wudx^j + C'(£)||V «||| 3 + 2e||V 2 i7||^ 2 . 

Similarly, we have 

H-Vu- Hdx^j + C(e)\\Vu \\ 3 L3 + 2 e\\V 2 H \\ 2 l2 . 
Putting (|3.9j) - (j3.13l) into (13.71) . and recalling (13.11) . (12.111) yields 
B\t) + J p\u\ 2 dx 

< C'||Vu || 2 2 +U(e)||Vn ||| 3 +4e||V 2 i7 || 2 2 

< C(s) (||Vu ||| 2 + l) + C'|||tf||Vff|||£ a +e\\p 1 / 2 u \\ 2 L2 +4 £ ||V 2 i7||| 2 . 

where 

B(t) =^||Vu ||| 2 + ^ y — ||divu || 2 2 - J divuPdx 
— ^ j divu\H\ 2 dx + J H • Vu • Hdx. 

Next, multiplying (|l.ll ) 2 by Ai7, and integrating by parts over M 2 , we have 


_d 

di 


\VH\ 2 dx + 2v / \V 2 H\ 2 dx 


<C j \Vu\\VH\ 2 dx + C J \X7u\\H\\AH\dx 

< C||Vu|| L 2||Vi7|| 2 4 +C||V'«|| L 2||i7|| L oo||V 2 /7|| i 2 

< C\\Vu\\ L 2 \\VH\\ L 2 \\V 2 H \\ L 2 + C\\Vu\\ L 2 \\V 2 H 

< ^||V 2 F|| 2 2 + CWVuWl, + C'||Vu||| 2 ||Vi7|| 2 2 . 


L 2 


Choosing C\ suitably large such that 

^I|Vh|| 2 2 + ||Vi7|| 2 2 — C'||P|| 2 2 

< m + Ci\\\ 7 H \\ 2 L 2 < C\\Wu\\l 2 + C\\VH \\ 2 l2 + C\\Pf L2 , 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 













adding (13.161) multiplied by C\ to (13.141) . and choosing e suitably small lead to 
(B(t) + Ci||ViL||| 2 ) / + IJ ( p\u \ 2 + vCi\\ 7 2 H\ 2 ) dx 
< C + C\\\H\\VH\\\ 2 l2 +C'(||V «||| 2 + \\vh\\ 4 l2 ). 

Integrating (13.181) over (0,T), choosing q = 4 in (|3.3D . and using (13.21) and Gronwall’s 
inequality, we obtain (13.51) . We can get (13.61) immediately by ( 12 . 111 ) and (13.51) . 

□ 


Next, we get some basic energy estimates on the magnetic field H. 
Lemma 3.3 Under the condition (13.11) . it holds that for 0 < T < T*, 

sup \\\H\\VH\\\ 2 l2 + I (\\A\H\ 2 \\ 2 L2 + \\\AH\\H\\\ 2 L2 )dt<C. 

Proof: We will follow an idea in [25]. For 01,02 £ {—1,0,1}, denote 
H(a±,a 2 ) = a\H l + o^iP 2 , tt(ai, 02 ) = a\u l + a 2 U 2 . 


It thus follows from (11.11) -? that 

H t - uAH = H ■ Vu - u • VF + Hdivu. 

Integrating (13.21[) multiplied by 4i/ _ 1 i7A |^| 2 over R 2 leads to 


-1 


VI H 


2112 

IIl 2 


+ 2||A|F| 2 1 | 2 2 


= 4 / \\/H\ 2 A\H\ 2 dx-du- 1 / H- X7u-HA\H\ 2 dx 


+ 4z/ 


-1 


■7 


u- V\H\ 2 A\H\ 2 dx 


d\\u\H\ 2 A\H\ 2 dx + 2 v 

< C'llVulll. + C'HViflll. + C|||i/| 2 ||l 4 + ||a |^| 2 || 2 2 

where we have used the following simple fact that 


7 


u- V\H\ 2 A\H\ 2 dx = -2 J diU-X7\H\ 2 di\H\ 2 dx + J divu\V\H\ 2 \ 2 dx 

<C\\Vu\\% + C\\\/H\\ 4 Li + C\\\H\ 2 \\ 4 L ,. 

Integrating (I3.22[) over (0,T), and using (13.31) and (13.61) . we obtain 


sup ||V|i^| 2 1|? 2 + 

0 <t<T 


\\A\H\ 2 \\ 2 L2 dt < C. 


(3.19) 


(3.20) 

(3.21) 


(3.22) 


(3.23) 


Noticing that 

\\\\7H\\H\\\ 2 l2 <||V|F(1,0 )| 2 ||| 2 + ||V|^(0,1 )| 2 ||| 2 

+ l|V|FT(l, 1)| 2 || 2 2 + ||V|iL(l, —1)| 2 ||| 2 , 

and that 

ll|A^||/f|||l 2 <C||V^||l 4 + ||A|^(l,0 )| 2 ||| 2 + ||A|#(0,l)| 2 ||i 2 
+ l|A|^(l, 1)| 2 || 2 2 + ||A|i7(l, —1)| 2 || 2 2 , 
then we can get (13.191) from (13.231) - (13.251) . 


(3.24) 


(3.25) 

□ 


In order to improve the regularity estimates on p , u and H, we start some basic 
energy estimates on the material derivatives of u. 
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Lemma 3.4 Under the condition (13.11) . it holds that for 0 < T < T*, 


sup (||p 1 / 2 u||i 2 + ||V< 4 )+/ ||V«||i a dt < C. 

<u+<ut V / ./n 


0 <t<T 


(3.26) 


Proof: We will follow an idea due to Hoff [9]. Operating d/dt + div(u • ) to (II. lD ^ 


and multiplying the resulting equation by v ?, one gets by some simple calculations that 
i p\v?\ 2 dx^ = ^J + d\w{uAu^))dx 

+ (p + A) / ii J (dtd-jidWu) + div(udj(divu)))dx 


— J v?(djP t + div(udjP))dx 
— v?(dtdj\H \ 2 + div(udj\H\ 2 ))dx 

5 

+ / u\d t {H ■ V&) + div(u(H ■ VW)))dx = ^Ji- 

J i=1 

First, integration by parts gives 

J\ = p J v? (Auj + div(uAu^))dx = —p J (d t u 3 dju\ + A nPu ■ Wv?)dx 
= —p [ (|Vh| 2 — did?u k dkdiui — d,u 3 dpi k dpu? + A vPu • X7u 3 )dx 


= —p I (|Vu | 2 + diU J dku k dpi 3 — dpi 3 diU k dkU J — diU 3 diU k dkU 3 )dx 


- +C , H Vn lll 4 - 


Similarly, 


J 2 < -^4^l|divft ||| 2 + C||Vu||^ 4 . 


(3.27) 


(3.28) 


(3.29) 


It follows from integration by parts, and (13.81) . (13.51) that 

J3 = — J u\djPt + di v(udjP))dx = J {djV?Pt + djPu ■ WvP)dx 
= — J ((7 — l)(9 J -u :? Pdivn + djii 3 div(Pu) + Pdj(u • X7u 3 ))dx 
= — J ((7 — l)djU 3 Pdivu + djii 3 dxv{Pu) + Pu ■ Vdjii 3 + PdjU • S7vP)dx (3.30) 
= — j ((7 — l)djiPPddvu + PdjU ■ Vu 3 )dx 
<^\\Vu\\l 2 +C. 
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Next, it follows from (I1.1IL and CD that 


J 4 = / djU 3 H • H t dx + - I u, ■ v u, wj 


djidd\vu\H\ 2 dx — — j L/j i 


u ■ Vvjdj\H\ 2 dx 

diU ■ '\7v?\H\ 2 dx 

J 11 2 J 

+ J djV?H • ( H • Vu + z^A H — Hdivu)dx 

<C j \Vu\\S7u\\H\ 2 dx + C J \Vu\\AH-H\dx 

< f l|VA||i 2 + C||Vn|| 4 i4 + C\\H\\ 8 l8 + C\\AH ■ H\\\ 2 . 

Similar to (13.311) . we estimate J 5 as follows 

J 5 < ^l|VA||| 2 + C||Vn|| 4 i4 + C\\H\\ 8 L s + C\\AH ■ H\\ 2 l2 . 


(3.31) 


(3.32) 


Putting (I3.28I) - (I3.32I) into ()3.27D . which together with (13.31) gives 

^11 P 1/2 n\\h + l|V«||i 2 dx < C||W|| 4 4 +C\\\AH\\H \\\ 2 l2 +C. (3.33) 

Integrating (|3.33l) over (0,T), using (13.61) and (13.191) . one has 

sup ||/9 1 / 2 ?1||^2 + [ ||VL|| 2 2 dt < C. (3.34) 

0 <t<T JO 

Then (13.26(1 can be obtained directly from (12.111) . (13.11) . (13.21) . (13.191) and (13.341) . 

□ 


Next, the following Lemma 13.51 combined with Lemma 12.41 will be useful to estimate 
the L p -norm of pu and obtain the regularity estimates on p. 


Lemma 3.5 Under the condition (133]), then there exists a positive constant N\ de¬ 
pending only on Nq, T and the initial data such that for 0 < T <T*, 



/ p(x,\ 

t)dx > j. 

(3.35) 


Ib Ni 

4 


Proof: First, multiplying (11. 1 D 1 

by (l + |x 

2 ) 4 / 2 and integrating the 

resulting equality 

over M 2 , we obtain after integration by parts and using both (13.21) and the fact that 

/ 

p dx = J 

po d.x = 1. 

(3.36) 

This leads to 




sup 

0 <t<T 

jpd+\ 

x\ 2 ) 1/2 dx < C. 

(3.37) 


Next, for N > 1, let <pn(x) be a C'q°(M 2 ) function such that 

0 < Vn(x) < 1, ¥>*(*) = { J !f|x|>^V, |V^|<2iV- 1 . 
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It follows from (II. 11^ . (13.21) and (|3.36D that 


d_ 

dt 


pp^dx = 


pu ■ Vp^dx > —2 N 1 



> - CN- 1 


which implies 


J pp^dx > 


pop^dx — CN l T. 


It thus follows from (11.111) and (13.381) that for N\ = 2(2 + Nq + 4 CT), 


(3.38) 


/ pdx > / 

J b Ni J 


pp^dx > 

2 4 


which gives (13.351) . The proof of Lemma 13.51 is completed. 


□ 


The next key lemma is used to bound the density gradient and L 1 (0, T; L°°)-norm 
of Vu. 


Lemma 3.6 Under the condition for any q > 2, it holds that 

SU P GlPlIiPrw 1 ^ + ||Vtt||^i) + f \\y 2 u\\f q dt < C. (3.39) 

0 <t<T JO 

Proof: In fact, for p € [2 ,q\, |Vp| p satisfies 

(\\7p\ p ) t + div (\\7 p\ p u) + (p — l)|Vp| p dmi 
+ p\Vp\ p ~ 2 (VpYVu(Vp) + pp\Vp\ p ~ 2 Vp ■ Vdivu = 0. 


Thus, 

|||Vp|| LP < C( 1 + ||Vu||l°o)||Vp||lp + c\\v 2 u\\ lp (3 . 40) 

< C(l + \\Vu\\ L ^)\\Vp\\ L P+C\\pu\\ LP + C\\\H\\VH\\\ LP , 

due to 

||V 2 u||lp < C (UphllLP + ||VP||lp + |||iL||VlL|||Lp), (3-41) 

which follows from the standard L p -estimate for the following elliptic system: 

pAu + (p + A)Vdivu = pu + VP + -V|iL | 2 — H ■ VH, u —>• 0 as |a:| —^ oo. 


Next, the Gagliardo-Nirenberg inequality, (12.81) and (12.91) implies 

||divu|| L oo + IMI L°° 
if C\\F\\ L 00 + C\\P\\ L 00 + C | | CO | | 

< C{q) + C7(g)||VP||lf (9 - 1)) + C(g)||Vu;||i/ (2(9 - 1)) 

< C(g) + C(q) {\\pu\\ Lq + \\\H\\VH\\\ Lq ) q /^-^ , 
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which, together with Lemma 1127771 (13.411) and (13.51) . yields that 

||Vu||l» (||divttH^oo + ||w||£,oo) log(e + ||V 2 'u||l9 ) + C 1 ]] Vit ||£2 + C 

<C (l + \\pu\\ q / q i2ig ~ 1]) + |||fl'||Vfl'|||^ (2(,_1)) ) 

• log(e + ||p«||l? + |||i?||Vi/|||L9 + 11 Vp 11 li ) 

<C (1 + ||pu|U« + ||| fr || Vfl -||| L ,) log(e + IIVpIIlO- 


It follows from Lemma [274j Lemma 13.51 (13.371) and ()3.26[) that, 



dt < C 



h + ||Vuf i2 ) < c. 


(3.44) 


Moreover, we have by Holder’s inequality, Gagliardo-Nirenberg inequality that, 


H\\VH\\\ Lq < C||fl-|| L 2,||V/f|| L 2, < C\\H\\ l ^\\VH\\ L 2\\y 2 H\\% 1)/q 
<C\\VH\\ l ,{\\H\\l2 + \\V 2 H\\ L 2), 


(3.45) 


integrating this over (0, T), which together with (|3.2D and f|3.5|> . yields 



\H\\VH\ 


2 

M 


dt < C. 


(3.46) 


Then, substituting (|3.43l) into (13.40[) where p = q, we deduce from Gronwall’s in¬ 
equality, (I3.44p and (I3.46P that 


sup \\Vp\\ L i < C, (3.47) 

0 <t<T 

which, along with (I3.43p . and (I3.4ip . shows 

(||V«|| L oo + \\V 2 u\\ 2 Lq ) dt < C. (3.48) 



Finally, taking p = 2 in (13.401) . one gets by using (13.31) . (13.51) . (13.481) . and Gronwall’s 
inequality that 

sup ||Vp|| L 2 < C. (3.49) 

0 <t<T 

The standard L 2 -estimate for the elliptic system, ()3. 191) . (|3.26l) and (|3.49l) lead to 

sup ||V 2 n|| 2 2 < C sup (||pft||| 2 + ||VP|| 2 2 + |||iL||Vtf||] 2 2 ) < C, (3.50) 

0 <t<T 0 <t<T 

which together with (13.361) . (13.11) . (13.51) and (j3.47p ~ fl3.49j) finishes the proof of Lemma 

[3761 □ 


Next, it follows from (12.31) . (13.351) . and the Poincare-type inequality |H Lemma 3.2] 
that for s > 2, p £ (0,1], and t € [0, T], 

||ux _1 || L 2 + Wux-^W^/r, < C(s,r])\\p 1/2 u \\ L 2 +C(s, 77)11 Vrt|| L 2 , (3.51) 

which together with (13.21) and (13.51) gives 

||ux 1 11 j ^2 + \\ux~ v \\ La /ri < C(rj,s). (3.52) 

With the help of (|3.52p . we can get the following spatial weighted mean estimate of the 
density, which has been proved in G2H Lemma 4.2]. 
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sup \\x a p\\ L i nH i nW i, q < C. 

0 <t<T 


Lemma 3.7 Under the condition (13.11) . it holds that for a > 1, q > 2 and 0 < T < T*, 

(3.53) 

Lemma 3.8 Under the condition dm it holds that for a > 1 and 0 < T <T*, 

r T 

sup \\Hx a/2 ||| 2 + / \\VHx a/2 \\ 2 L2 dt < C, (3.54) 


0 <t<T 


r T 

sup \\VHx a/2 \\ 2 L2 + / \\V 2 Hx a/2 \\ 2 L2 dt < C. 

0 <t<T Jo 


(3.55) 


Proof: First, multiplying (II.ID s by Hx a and integrating by parts yields 


| H\ z x a dx ) +u \X7H\ z x a dx = - / \H\ z Ax a dx 


+ / H ■ Vu ■ Hx a dx - 


c a dx - i f dWu\H\ 2 x a dx + - / \H\ 2 u ■ Vx a dx = J2 K ^ 
•' J i =i 

Direct calculations and (13.51) yield that 

|iFi| <cj \H\ 2 x a x ~ 2 log 4 (e + \x\ 2 )dx <C j \H\ 2 x a dx, 


(3.56) 


(3.57) 


and that 

|iF 2 | + \K 3 \ <C j\Vu\\H\ 2 x a dx < C\\Vu\\ L2 \\Hx a/2 \\ 2 L4 

< C\\Vu\\ L 2 \\Hx a / 2 \\ L 2 (\\VHx a / 2 \\ L 2 + \\HVx a / 2 \\ L 2 ) 

< C\\Hx a/2 \\ 2 L2 + j\\VHx a / 2 \\ 2 L2 . 

Then, it follows from Holder’s inequality, (12.11) and (13.521) that 

\K 4 \ < C\\Hx a / 2 \\ L 4Hx a / 2 \\ L 2\\ux- 3 ^\\ L 4 

< C\\Hx a / 2 \\ 2 L 4 + C\\Hx a / 2 \\ 2 L2 ||ux _3 / 4 ||^ 4 

< C\\Hx a/2 \\ 2 L2 + j\\VHx a / 2 \\ 2 L2 . 

Putting (13.571) - (13.591) into (13.561) . after using Gronwall’s inequality, we have 


(3.58) 


(3.59) 




sup x a \H\ z dx+ I x a \VH\ z dxdt < C. 
0 <t<T J Jo 


(3.60) 


Next, multiplying (II.IIP by A Hx a , integrating the resultant equation by parts over 
2 , it follows from the similar arguments as (|3.16l) that 


| VH\ z x a dx) +u j \AH\ z x a dx 


<C J \VH\\H\\S7u\\Vx a \dx + C J \VH\ 2 \u\\Vx a \dx + C J \VH\\AH\\Vx a \dx 

5 


+ 


C f \H\\X7u\\AH\x a dx + C f \Vu\\X7H\ 2 x a dx = L t . 

i =1 


(3.61) 
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Using Gagliardo-Nirenberg inequality, (13.261) and (13.541) . it holds that 


Li <C j |Vi/'||fr||V«|x a (x _1 |Vx|)ds 

<C\\Hx a / 2 \\ A L , + C||V«|| 4 4 + C\\'VHx a ^ 2 \\ 2 L 2 

<C\\Hx a/2 \\ 2 L 2 (j|ViLx a/2 || 2 2 + \\Hx a/2 \\ 2 L 2 ^j + C + C\\VHx a / 2 \\ 2 L 2 

<C + C\\VHx a / 2 \\ 2 L2 , 


L 2 <C J |Vi/| (4a " 1)/(2a) x (4a " 1)/4 |Vi7| 1/(2a Vl®“ 1/2 S _1/4 |Vx|dx 

< C|||Vi?| (4a - 1)/(2a) x (4a - 1)/4 || i3 |^|||V/f| 1 / (2a) || L 8 a ||^" 1 / 2 || L8a ^ 3 63 ^ 

< C\\X7Hx a/2 \\ 2 L2 +C\\VH\\ 2 l4 

< C\\VHx a / 2 \\ 2 L2 + C||V 2 tf || 2 2 + C, 

L3 + L4 < ||A^x a/2 || i 2||V//x a/2 || L 2 + ||A^x“/ 2 || L 2||iJx°/ 2 || L 4||Vw|| L 4 

< ^||Aiix a/2 || 2 2 + C\\VHx a / 2 \\ 2 L2 + C, ^' 64 ^ 

L 5 <C||V'u|| L co||V//x a/2 || 2 2 

<C||Vw||^V 2)/(29 “ 2) ||V 2 u||^ (29_2) ||Vi?x a/2 || 2 2 (3.65) 

<C{l + \\V 2 u\\ 2 Lq )\\VHx a / 2 \\ 2 L2 . 

Noticing the fact that 


|V 2 iL| 2 x a dx = J \AH\ 2 x a dx — J didkH ■ dkHdiX a dx 

+ J didiH ■ dkHdkX a dx 

< f \AH\ 2 x a dx + l [ \X7 2 H\ 2 x a dx + C f \VH\ 2 x a dx. 


(3.66) 


Submitting (I3.62[) ~ (13.65D into (13.611) . and using (|3.66|> . we obtain 


iViLlVdx 


+ / |V 2 iL| 2 x a dx 


< C( 1 + ||V 2 «||| 9 )||Vi/x a / 2 || 2 2 + C(||V 2 tf || 2 2 + 1), 


(3.67) 


which together with Gronwall’s inequality, (13.51) and (13.391) yields (13.551) . The proof of 
Lemma 13.81 is finished. □ 


Lemma 3.9 Under the condition (13.11) . it holds that for 0 < T <T*, 
sup (\\p 1/2 u t \\ 2 L2 + ||#f || 2 2 + l|Vi7||^i) 

0 <t<T v ' 

+ [ T (l|V ^|| 2 2 + ||Vtf t ||| 2 + ||v 2 i7||| 9 ) dt < c. 
Jo 


(3.68) 
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Proof: First, the combination of (13.531) with (13.521) gives that for any p E (0,1] and 
any s > 2 , 

\\p v u\\ L s/ v + Wux^Whs/r, < C(r), s ). (3.69) 

Differentiating (II.lib with respect to t gives 


pu tt + pu ■ Vu t — pAu t - {p + A)Vdi vu t 
= — pt(pit + u ■ Vu) — put ■ Vu — VP t + (h • VH — -V\H\ 2 


(3.70) 


Multiplying (13.701) by u t , then integrating over R 2 , we obtain after using (I1.1D 1 that 

J P\ Ut \ 2dx + J (HVut| 2 + (p + A)(divu t ) 2 ) dx 

= —2 / pu ■ Wut ■ utdx — / pu • V(u • Vu • ut)dx 


— j pu t ■ V u ■ u t dx + j Pfdivutdx + / ( H ■ VH — -V|iL | 2 ) utdx 


(3.71) 




i= 1 


Similar to the proof of }20l Lemma 4.3], and using (13.391) . we have 
Jj < ^\\^ u t\\ 2 L 2 + C ( ll^^llzp + ||p 1/2 u t || 2 2 + 1 

2—1 

< ^l|Vut||| 2 + C (||p 1//2 Ui|| 2 2 + l) • 

For the term J 5 , we obtain after integration by parts that 


(3.72) 


J 5 = - / H, ■ V», ■ Hdx - j H - V Ui H t dx + J H ■ H,d iv U(d x 

<C||/7|| L 4||^|| L 4||Vu t || i2 

< ^l|V«t||ia + <J||Vflt|||,2 + C(5)||F7 i || 2 2. 


(3.73) 


Substituting (13.721) and (J3.73D into (13.711) leads to 

^\\p 1/2 ut\\b + l|V^||| 2 < C 5\\V Ht \\ 2 L 2 + C(6) (||p 1/2 u t ||i 2 + \\H t \\ 2 L2 + 1 ) • (3.74) 

Next, differentiating with respect to t shows 

H tt — H t ■ Vu — H ■ Vut + ut • V H + u • V H t + H t divu + Hdivut = vAH t . (3.75) 

Multiplying (|3.75l) by H t , then integrating the resulting equation over R 2 , and using 
(13.31) . (13.51) . (13.511) and (13.541) . we have 

1 d 


2 dt 


•J \VH t \ 2 dx = J 


H t fdx + 1 / / \VH t rdx = I H t ■ Vu • H t dx 


- 1 / divulff, \ 2 dx + J H-Vu, . iMx + f u, ■ VH, ■ Hd X 

<C||^|| 2 4 ||Vu|| i2 +C7||77|| L 4||iL t || L 4||Vu t || L2 

... , _1_|. 11 / a . 2a —1 _ 2a —1 

+ ( 7 || |-ff| 2a ||^8a || (H X 2 ) 2“ || 4a ||UfX 4 11 j^Sa \ \ Ht \ \ ]? 

L ‘^ a ~ i 

< ^l|VL7t|| 2 2 + C(5) (||^||1 2 + ||p 1/2 ^||i 2 + ||Vu t ||| 2 ) , 


(3.76) 
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which gives 
d 


-^WHtWv + ||Vtf t ||| 2 < C\\H t \\ 2 L 2 + C\\p 1/2 u t \\ 2 L 2 + C 2 \\Vut\\b- (3-77) 


Then adding (|3.74|) multiplied by C 2 + 2 to (13.771) , choosing <5 suitable small, we obtain 



(3.78) 


which together with Gronwall’s inequality yields 



(3.79) 


sup 

0 <t<T 


Finally, it follows from (13.791) . (13.31) . (|3.52p . (I3.26p . (13.551) and Gagliardo-Nirenberg 
inequality that 


V 2 H\\ 2 L 2 < C\\\u\\VH\\\ 2 l2 + C\\\H\\Vu\\\ 2 l2 + C\\H t \\ 2 L 2 

< C\\ux~ a ^\\ 2 L s, HI V77| 1 / 2 x a / 4 ||| 4 1| |VT/) 1 / 2 1|| 8 
+ C\\H\\ 2 l 4Vu\\ 2 l4 + C 

< C\\VHx a/2 \\ L 2\\VH\\ Li + C 


(3.80) 



and 



_|_ II TJ 11^/9 
+ \\ H t\\ L 2 


VH t \\ 2 ^~ 2)/q )dt 



+ UVull^i + ||ViZt||£ a + l)dt 



(3.81) 


which combined with (13.801) leads to (13.681) . Thus the proof of Lemma 13.91 is finished. 

□ 

Now we are in a position to prove Theorem 11,11 

Proof of Theorem 11.11 Suppose that (11.141) were false, that is, (13.11) holds. Note 
that the generic constant C in Lemma 13.21 and Lemma 13.6113.91 remains uniformly 
bounded for all T < T*, so the functions (p,u,H) = lim^y* (/?, u, H)(x, t) satisfy 
the conditions imposed on the initial data (11.121) at the time t = T*. Furthermore, 
standard arguments yield that pu G C([0, T]; L 2 ), which implies 


(pu)(x,T*) = lim (pu) G L 2 . 
t—*T* x 
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Hence 


—fj.Au — (fj + A)Vdivu + VP — (V x H) x Pd\ t=T * = \fp(x, T*)g(x), 


with 


g(x) 


p(x, T*)-^ 2 (pu){x, T*), for x € {x\p{x, T*) > 0}, 
0, for x € {x\p(x, T*) = 0}, 


satisfying g £ L 2 due to (|3.26l) . Thus, ( p,u)(x,T*) satisfies (| 1.13 f) also. Therefore, one 
can take (p,u, H){x,T*) as the initial data and apply Lemma 12.11 to extend the local 
strong solution beyond T*. This contradicts the assumption on T*. We thus finish the 
proof of Theorem 11.11 

Acknowledgments. The author would like to express great gratitude to Prof. 
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